Breathers in a system with helicity and dipole interaction. 
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Recent papers that have studied variants of the Peyrard-Bishop model for DNA, have taken into 
account the long range interaction due to the dipole moments of the hydrogen bonds between base 
pairs. In these models the helicity of the double strand is not considered. In this particular paper 
we have performed an analysis of the influence of the helicity on the properties of static and moving 
breathers in a Klein-Gordon chain with dipole-dipole interaction. It has been found that the helicity 
enlarges the range of existence and stability of static breathers, although this effect is small for a 
typical helical structure of DNA. However the effect of the orientation of the dipole moments is 
considerably higher with transcendental consequences for the existence of mobile breathers. 
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I. INTRODUCTION 

A great deal of attention has been paid to the interplay 
between geometry and nonlinearity in locating problems 
in recent years. The relationship between geometry and 
nonlinearity has an important role in the functions of 
some biomolecules, such as DNA, where the localization 
of energy has been put forward as a precursory mecha- 
nism of the transcription bubble Q, and moving local- 
ized excitations as a method of transporting information 
along the double strand Q . 

The fact that hydrogen bonds that link each pair of 
bases in DNA have a finite dipole moment, has brought 
about the introduction of models [^-^| with long range 
dipole-dipole interaction. Apart from its theoretical in- 
terest, this interaction becomes relevant when the sec- 
ondary structure of DNA is considered. The shape of 
the molecule can influence the localization and transport 
properties of energy, which is thought to play a biological 
function Q). Some of these models || study the effects 
of the curvature in a chain of nonlinear oscillators us- 
ing the discrete nonlinear Schrodingcr equation. Other 
models consider Klein-Gordon systems to study kinks j|, 
breathers in curved chains j^] or breathers with two com- 
peting interactions Pj. However, all these models with 
long range interaction fail to take into account the pe- 
culiar helicoidal structure of the DNA chain, although 
this has been considered in some models H without the 
dipole interaction. 

In this paper, we study the effect of helicity on the 
properties of breathers in a Klein-Gordon model with 
dipole-dipole interaction. These periodic nonlinear lo- 
calized oscillations in discrete systems are very localized 
excitations that appear as a consequence of the nonlin- 
earity and discreteness of the system g). They are spe- 
cially suitable for biomolecules when considering excita- 
tions that involve a few units, that is, far from the con- 
tinuous limit. They can be static but, under certain con- 



ditions, also move and transport energy along the system 

We have found that the introduction of helicity en- 
hances the stability of static breathers, although this ef- 
fect is relatively small for the typical helicoidal structure 
of the DNA. On the other hand, the profile of the static 
breathers and the properties of moving ones are strongly 
dependent on the relative orientation between the dipole 
moments. 



II. THE MODEL 



The model is inspired by the primary structure of 
DNA, with dipole moments perpendicular to the helix 
axis, and where the stretching of the hydrogen bonds 
within base pairs is described as a variation of the dipole 
moments. More detailed justification of the model can 
be found in 101 . 




FIG. 1. Sketch of the model at equilibrium. The arrows 
represent the dipoles moments, perpendicular to the helix 
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We denote <p n the angle of the n-dipole with respect to 
a reference axis perpendicular to the helix axis. Then, 
the angle between the nearest neighbouring dipoles is 
Otw = 4>n — fin— I- We have considered this neighbouring 
angle constant along the chain, and it will be called the 
twisting angle. Thus, </>„+,„ — <f> n — n%6t w , and, therefore, 
27r /Otw dipoles arc needed to complete a turn of screw. In 
DNA, for example, the twisting angle is 36° and a turn of 
screw requires ten base pairs. Figure [l] shows a sketch of 
the model, where it can be appreciated that the system 
of dipoles have an hclicoidal structure. 

In the appropriate dimensionless variables, the Hamil- 
tonian of our system becomes 

N 



n+N/2 



^n^m 
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co$[9 tw (n - rn)] ) , (1) 



where N is the number of variables. The variables 



A" 



whe 



re u n±N 



represent, in the context 



of the Peyrard-Bishop model for DNA |L| , the transver- 
sal displacements of the two complementary nucleotides 
in the n-th pair with respect to the molecular axis. In 
our model, they describes the stretching of the dipoles 
with respect to their equilibrium length. V(u n ) is the on 
site potential, which, in DNA models, describes the hy- 
drogen bonds linking the two bases, and the parameter 
J measures the strength of the long range dipole-dipole 
interaction. We have chosen the on-site potential as the 
Morse potential, given by 



V{u n ) = \{e—-lf 



(2) 



The reason for this, is that it is a suitable potential 
for representing chemical bonds, being asymmetric, with 
a hard part, modeling the repulsion between atoms or 
molecules, and a soft part that becomes flat, modeling 
the breakage of the bond. 

The dynamical equations become 



Un + V'(u n ) + J 



n+N/2 

cos[O tw (n -m)\ 

' In — m| 3 

m=n-JV/2 1 1 



0, (3) 



where n=l ... N. To study the linear modes of the sys- 
tem we replace V' (u n ) in equations |^ with the linear term 
U n , which implies that the time has been scaled so that 
the linear frequency u)q = 1. Considering solutions of the 
form u n = e l 9™-™* the following dispersion relation is 
obtained: 



w k = 



\ 



N/2 



1 + 2JJ2 



cos(m9 tw ) 



cos(mq k ) (4) 



where = , with k = 1 . . . N due to the periodic 
boundary conditions. 



The variation of the phonon band with the helicity 
is shown in Figure ^, where the frequencies of the lin- 
ear modes are represented as a function of the twisting 
angle, 9 tw , for a fixed value of the coupling parameter 
J = 0.1. The effect of the twisting is a narrowing of the 
phonon band, which will enhance the range of existence 
and stability of the breathers. This has been confirmed 
numerically. 
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FIG. 2. Effect of helicity on the breadth of the linear spec- 
trum. The curves represent the lower and higher limits of the 
phonon band as a function of the twisting angle in degrees 
for fixed coupling parameter J = 0.1. to is in dimensionless 
units. 



III. BREATHER EXISTENCE AND STABILITY 



We have studied the existence and stability of 
breathers in this model using the standard numerical 
methods described in Ref. JTl| . 

The Morse potential is a soft potential with the conse- 
quence that the frequency of a breather has to be lower 
than the linear frequency too — 1- Thus, we have chosen, 
u>b = 0.8 so that the nonlinear effect will be significant 
but on the other hand not overtly strong, as the nonlin- 
earities in DNA are thought to be weak. 

First of all, the helicity influence the breather pro- 
file. As is shown in figure ||, for a fixed value of the 
coupling parameter J, the increase of the twisting angle 
produces a transition from a zigzag profile (the nearest 
neighbour oscillating in antiphase) to a bell profile (all 
dipoles oscillating in phase). This effect follows from the 
spatial profile of the phonon state with the lowest fre- 
quency since the breather frequency is below the phonon 
band, and all the higher harmonics are way too high to 
be relevant. For tw < ir/2 the interaction is effectively 
"antiferromagnetic" which leads to staggered phonons at 
the lower band edges. In the same way for 9 tw > it/2 
a "ferromagnetic" interaction is present which leads to 
a nonstaggered phonon at the lower band edge. The 
breather bifurcates from the lower band edge phonons 
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and thus retains the property of the phonon structure. 
For tw = it/2 the system separates into two noninter- 
acting sublattices: even and odd sites. As a result in this 
case, the nearest neighbors are at rest and the odd site 
sublattice remains uncxcited. 
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FIG. 3. Profiles of the one-site breather when twisting is 
increased for fixed coupling J = 0.1 and Morse potential. Di- 
mensionless units. 

One-site breathers are stable at low coupling as was 
proved by Aubry [Q. For any value of the twisting an- 
gle 9 tw < 90° they can be continued from the anticon- 
tinuous limit till Wb enters the phonon band. Just before 
the breather disappears, it becomes unstable due to the 
occurrence of a harmonic bifurcation in the evolution of 
the Floquet eigenvalues. The increase of the twisting en- 
hances the stability as is shown in Fig. [| (circles). This 
can be understood if we consider only the nearest neigh- 
bor interaction (NNI). Then the influence of helicity on 
the stability of the breathers could be described by an ef- 
fective coupling J e f f = J cos 6tw ■ The one-site breather 
without twisting lose its stability for a coupling value of 
J®. With twisting and only NNI this would occur for 
J c = Jc/cos8 tw (dash lines in Fig. ||), which concurs 
with the numerical results. 

The two-site breather, which consist of two neighbor- 
ing oscillators excited in phase, is also stable at low cou- 
pling. This can be understood in terms of Aubry's band 
theory . When coupling is increased a bubble of insta- 
bility appears due to Krein crunches between the phonon 
band eigenvalues and a localized eigenvalue of the Flo- 
quet operator. If we continue increasing the coupling the 
double breather definitely becomes unstable due to the 
occurrence of a subharmonic bifurcation. Again, the ef- 
fect of the twisting is to enlarge the range of stability 
toward higher values of the coupling parameter (full cir- 
cles in Fig. ||). This suggests that twisting might be a way 
to control the stability of the breathers in real systems. 

We have not considered the two-site breather in an- 
tiphase because it coincides with the one-site breather 
with zigzag profile, i.e., the Newton method converges 



to the same solution if we start at the anticontinuous 
limit with one non linear oscillator or with two nearest 
neighbor oscillators in antiphase. 
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FIG. 4. Range of stability of the one-site (circles) and two- 
sites (full circles) breathers as a function of the twisting angle 
in degrees. J c is the maximum value of the coupling parame- 
ter for which the breather is stable. The dash lines represent 
the values calculate within the NNI approximation. J c is in 
dimensionless units. 

A rather different situation is the one with tw > 90°. 
First, the one-site breather is always stable until it dis- 
appears. Second, the two-site breather is unstable at 
low coupling but becomes stable just before its extinc- 
tion. This behavior has important consequences for the 
mobility of these breathers as shown in the next section. 

For the sake of thoroughness we have also studied 
the effect of twisting with a hard </> 4 potential V(u n ) — 
v% + and a breather frequency iv^ = 1-2. Qual- 

itatively the results are similar except for the fact that 
breathers with 9 tw < 90° and breathers with 6 tw > 90° 
exchange their properties. 



IV. MOBILE BREATHERS 

Static breathers under certain conditions can be 
moved. The standard method to move a breather con- 
sists in perturbing its velocities with an spatially anti- 
symmetric vector, called the marginal mode Typi- 
cally, this method works within a certain range of param- 
eters near an exchange stability bifurcation. This occurs 
when a one-site breather becomes unstable and a two- 
site breather does the opposite at a nearby point. 

We have looked for mobile breathers in our system 
both with a hard 4> 4 potential and with a Morse poten- 
tial, but we have only had success with Morse potential 
and "ferromagnetic" interaction, i.e., 6t w > 90°. In this 
particular case, we found a similar situation to a stability 
exchange and we were able to move the breather perturb- 
ing it with the unstable localized mode of the two-site 
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breather. This is an interesting result because this con- 
figuration is equivalent to a chain of antiparallel dipoles 
twisted 7r — 8 tw < tt/2. In fact, we can only expect par- 
allel dipoles in synthetic DNA. 

A useful concept for describing the breather movement 
is its effective mass. If the norm of the perturbation ve- 
locity is A, the kinetic energy added to the breather by 
the perturbation is E — A 2 /2. The resulting translational 
velocity of the breather, v, is found to be proportional to 
A p0[ . Thus, moving breathers can be considered as a 
quasi-particle with a mass of m*, which can be defined 
through the relation m*v 2 /2 = A 2 /2. 
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FIG. 5. Dependence of the effective mass of the mobile 
breathers with the coupling parameter for 6 tw = 180°. The 
blank squares correspond to mobile breathers obtained from 
static one-site breathers. Full squares are obtained from 
static two-site breathers. 

We have studied the dependence of the effective mass, 
m* , with the coupling J. Figure ^ shows the result for 
antiparallel dipoles (8t w — 180°). Two different behav- 
iors were obtained depending on the initial conditions. 
If we perturb the two-site breather, we observe that its 
effective mass increases monotonically with the coupling 
(full squares in Figure ||). This reflects the fact that 
the two-site breather becomes stable with increasing cou- 
pling. But if a static one-site breather is chosen as the 
initial configuration, a minimum of m* appears showing 
the existence of an optimal value of the coupling to move 
this breather (see blank squares in Figure g). We think 
that this minimum expresses a balance between the two 
opposite effects produced by an increase of the coupling 
on the stability of the one-site and two-site breathers. 
Similar results are obtained for other values 6t w > tt/2. 

V. CONCLUSIONS 

We have considered a system of oscillating dipoles with 
helicoidal structure in order to study the effect of helic- 
ity on the existence and properties of breathers. This 
study is motivated by the helicoidal structure of DNA, 



and the fact that it can be described by a reduced dy- 
namic where the only degrees of freedom are the stretch- 
ings of the hydrogen bonds between base pairs, which 
have a finite dipole moment. In our model, the helicity 
produces a narrowing of the phonon band, and an en- 
largement of the range of existence and stability of the 
breathers, although this effect is small for a typical heli- 
coidal structure of DNA. 

The effect of the orientation of the dipole moments, 
i.e., if the twisting angle is greater or not than 90 de- 
grees, is however considerably higher. In particular, we 
have only found mobile breathers with a Morse potential 
and 9 tw > it/2. Understanding the necessary conditions 
to move a breather is still an opened question today. 
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